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We consider tunneling of vortices across a superconducting film that is both narrow 
and short (and connected to bulk superconducting leads at the ends). We find that 
in the superconducting state the resistance, at low values of the temperature (T) 
and current, does not follow the power-law dependence on T characteristic of longer 
samples but is exponential in 1/T. The coefficient of 1/T in the exponent depends 
on the length or, equivalently, the total normal-state resistance of the sample. These 
conclusions persist in the one-dimensional limit, which is similar to the problem of 
quantum phase slips in an ultra-narrow short wire. 



I. INTRODUCTION 

Vortex tunneling across narrow thin films is a quantum effect that limits superconduc- 
tivity in these systems. It does not rely on any vortices preexisting in the sample (due, e.g., 
to temperature or a magnetic field) but occurs as a result of quantum fluctuations — either 
a vortex entering the sample from the outside or a virtual vortex-antivortex pair created 
in the sample. A supercurrent (if present) exerts a force on the vortex and, for certain 
types of vortex motion, the total work done by this force is nonzero. That means that an 
amount of energy is taken away from the supercurrent, i.e., there is dissipation — an electrical 
resistance. 

Frequently, one considers vortex tunneling at a strong enough current, so that the vortex 
can nucleate close to the boundary: the closer it nucleates, the shorter is the distance it has 
to tunnel, and the larger is the tunneling rate.-*^ However, for very narrow samples (starting 
perhaps with a few tens of nm), the rate may remain substantial even at weak currents, 
when the vortex tunnels the entire width. This is the case we consider here, motivated in 
part by its similarity to quantum phase slips in a genuinely one-dimensional (Id) wire, a 
topic of much current experimental research.— ^'^'^^^^^^^^^^^^^ 

While the core of a vortex can for our purposes be considered point-like, it is important to 
take into account the long-range disturbance that vortex motion produces. Low-dimensional 
superconductors support a gapless plasmon mode,— which in the presence of a nearby 
ground plate (the case considered here) has an acoustic dispersion law, with a speed Cq. 
This leads to a new length scale, Ip = Ticq/T (T is the temperature), to which other length 
scales can be compared. In particular, one can distinguish between long wires, those of 
length L ^ /p, for which the final state of plasmons does not depend on the boundary 
conditions at the ends, and short ones, L <C /p, for which it may. 

For long Id wires, plasmon production has been shown theoretically to significantly affect 
the tunneling raXeM^i^ In the presence of disorder, it leads to a power-law current-voltage 
dependence in the superconducting state at low temperatures,— V{I) oc J° (a state is 
superconducting if a > 1). Here, we consider the opposite, short- wire, limit.* This choice is 
motivated by a puzzle in the existing experimental results: while in a long wire the observed 



* We will often use the system of units with 7i = cq = 1, in which the short- wire condition is simply LT <^ 1. 
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nonlinear V{I) curve is indeed well described by a power law,— in short wires, no power 
law in current or in temperature has been detected.— Our results provide an explanation 
for that. 

A general theory that allows one to compute the effect of plasmons on vortex tunneling 
for samples of various sizes can be constructed along the following lines. As well known (and 
reviewed, for example, in Ref. ITol ). two-dimensional superfiuids have a dual description, in 
which vortices are viewed as charges and plasmons as "photons", so that the theory maps 
onto planar quantum electrodynamics (QED). We present a derivation of this, based on a 
path- integral identity, in Sect. [Tll 

The version of QED that we use in this paper is "quenched", in the sense that it con- 
siders only a single vortex and neglects interaction with additional vortices that may tunnel 
nearby. We expect this to be a good approximation as long as one stays away from a 
superconductor-insulator transition. (For Id wires, a mean-field-type theory that takes into 
account interactions between quantum phase slips has been recently proposed in Ref. [20I .) 
In the quenched limit, the computation of the plasmon action amounts essentially to a 
Euclidean (imaginary time) version of the classical radiation theory. (The imaginary time 
appears since we are considering a tunneling process.) An important aspect of the theory is 
formulation of the boundary conditions at the ends of the sample. This is described in Sect. 
Illlt and the solution for the plasmon field is given in Sect. HVl 

Our final result is that, in the superconducting state of narrow short wires connected to 
bulk superconducting leads, the resistance due to vortex tunneling, at small temperatures 
and currents, is no longer a power law but an exponential in 1/T. We trace this stronger 
suppression to a large gradient energy that the system must have already when it enters the 
classically forbidden region. The coefficient of 1/T in the exponent depends inversely on the 
length L. 

We wish to reiterate that this result applies only in the superconducting state, where the 
tunneling events (instantons) are rare, and does not preclude the possibility of a transition 
to an insulating state at larger tunneling rates. 

The tunneling process we consider here is in addition to and competes with the classical, 
over-barrier process. The latter is an analog of a thermally activated phase slip in the 
Id case.— As we will see, despite the above-mentioned suppression, the rate of vortex 
tunneling is exponentially larger than the rate of thermal activation over a broad range of 
parameters. We hope that the difference in both the magnitudes of the resistance and its 
dependence on the length will allow one to distinguish between the two effects experimentally. 

We should also note the difference in the starting points for the theories of these two 
effects. The energy of the LAMH saddle point^i^^ is due to the phase slip core (and conse- 
quently depends strongly on the Ginzburg-Landau coherence length ^). In contrast, in our 
case, the activation energy is that of the initial tunneling state, which lies far from the top 
of the potential barrier. Depletion of the order parameter in this state is still small, and we 
can use the phase-only theory, in which vortex cores are essentially point-like. 



The Lagrangian density of the phase-only theory that describes a superconducting film 
in the presence of a nearby ground plate is 



II. DUALITY MAP 




(1) 
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In the second equality here we have switched to the notation of special relativity, by using 
instead of time t the coordinate x° = Cot, where cq = yJgKg is the plasmon speed. (In 
addition, = x and x^ = U-) In what follows, we choose units of length and time so that 
Co = 1. Greek indices take values 0, 1, 2, and summation over a repeated index is implied. 

The field 9 is the phase of the order parameter, but no assumption is made about the 
existence of long-range order, i.e., we do no require the expectation value of exp(i^) to be 
nonzero. All that is required for superfluidity is that the stiffness Ks renormalizes to a 
nonzero value in the infrared. 

If 9 were a single- valued smooth function of (x, y), plasmons would be the only excitations 
in the system, and the theory would be completely Gaussian. To describe vortices, we allow 
9 to be multivalued. Alternatively, we could make it discontinuous by drawing explicit 
branch cuts at vortex positions, but we will be using the first approach. Then, V6' is smooth 
outside vortex cores. 

The phase-only description ([T]) does not resolve vortex cores, so a short-scale cutoff of 
order of the Ginzburg-Landau coherence length is implied. 

Let us remark on the issue of gauge invariance and the apparent absence of electromag- 
netic fields from ([1]). In thin films, magnetic field of a vortex extends over a large area, 
determined by the transverse screening length If (which in thin films can exceed 

100 /im) is much larger than the smallest dimension of the film, the magnetic field can be 
neglected, and it is possible to choose a gauge such that the vector potential A is close to 
zero. The remaining gauge freedom can be used to make the scalar potential Aq go to zero 
away from the film. This fixes the gauge completely, and if 9 denotes the phase in this 
gauge it is in effect gauge- invariant. Integrating out Aq produces a capacitive term, which 
eventually becomes the first term in ([T]). 

We have not included in ([1]) a "topological" term, proportional to dt9, that gives rise to 
the Magnus force on a vortex. We consider films that have significant amounts of disorder, 
and in disordered superconductors the Magnus force is small.— 

Eq. ([1]) allows one to describe dissipation of energy of vortex motion into plasmons, 
but does not include any dissipative mechanisms related to normal electrons at the vortex 
cores. This is justified in the limit of strong disorder, since transfer of energy to the normal 
component in this case is inhibited by the short electron mean-free path. 

Because in the presence of vortices 9 is multivalued, the expression 

where e^'^^ is the unit antisymmetric tensor, is nonzero and indeed is the vortex current ( J° 
is the vortex density). Setting 

qx = dx9 , (3) 

we can write the current ([2]) as 

J/. ^ ^e'^'^^d.qx . (4) 

The equation of motion following from ([T]) is 

d^dx9 = d^qx = . (5) 

[In the presence of a topological term, the temporal component of ([3]) is replaced by go = 
do9 + const., and eq. f|TT]) below is modified accordingly. The expression (jl]) for the current 
and the equation of motion ([5]) are both unaffected.] 
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In what follows we restrict our attention to configurations satisfying the equation of 
motion ([5]). In the real-time version of the theory, they describe motion of an arbitrary 
number of real vortices in the presence of supercurrents and plasmon waves. In the Euclidean 
(imaginary time) version, to which we turn shortly, solutions to the equations of motion will 
determine the most probable tunneling paths (instantons) responsible for the quantum decay 
of supercurrents. 

The real-time action corresponding to the Lagrangian density ([T]) is 

S=\K,j d'^ed^ed^x = J q^q^d'x , (6) 

where d^x = dx^dx^dx^. If is a solution of the equation of motion ([5]), there is a path- 
integral identity — the duality map: 

e'' = J Vf,,VX exp J iK, [-\r f,u - \t^'''f,.qp + \\t^'''d,f,^^ d'x , (7) 
where f^u is antisymmetric in fi, v and is subject to the boundary condition 

e^'"'f,.n,\^ = -2n% ; (8) 

Up is the normal to the boundary (6) of the spacetime volume. 

Because the path integral in ([7]) is Gaussian, the map can be verified directly. First, note 
that the path integral over A enforces the Bianchi identity 

e'^'dpf,. = . (9) 

Since in ([7]) ff^i, is an independent variable and not (yet) a curl of some gauge field, Qj is 
not really an identity but an independent equation of motion; however, we keep the familiar 
term. 

Next, integrating over /^,^ amounts to solving the saddle-point equation 

r + ^^"'{qp + dpX) = . (10) 

Taking curl of this and using eqs. and we obtain d'^dpX = 0. On the other hand, 
applying the boundary condition ([8]) in eq. ffTOl) tells us that the normal derivative of A at the 
boundary is zero. For a spacetime volume of a simple shape, these conditions are sufficient 
to reduce A to a constant, which then drops out of ffTOl) . Eq. ffTOl) becomes 

= -e^'^^pq^ = -e^^'^PdpO . (11) 

Substituting this back in ([7]) we confirm that the original action ([H]) is recovered. In terms 
of the saddle-point value flTTl) . this action can be written as 

S=\Ksj rfpud^x. (12) 

Now, let us take an alternative ("dual") view and solve the identity ([9]) explicitly by 
introducing a new gauge field (quite distinct from the electromagnetic potential A^): 



f pu dp,Ojv dyOj^ 



(13) 
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Differentiating (fTTI) . we obtain Maxwell equations for this field: 

dj^" = -27r , (14) 

where J'^ is the vortex current (jll). The usefulness of this dual view is that it allows one to 
determine /^j^, and the corresponding action, for any prescribed motion of vortices. Using 
( ITTl) . one can then find the derivatives of the phase 6. 

The dual theory also allows one to take into account the backreaction of produced "pho- 
tons" on the vortex motion, by extremizing the total action with respect to that motion itself. 
In the problem to which we apply this theory here, backreaction will determine the region 
where the vortex prefers to tunnel. (This region turns out to be the middle of the wire — in 
contrast to the resistively shunted Id case, where quantum phase slips occur preferentially 
near the ends.~) 

To any solution of eq. ( |T4l) we can add a homogeneous solution — a static uniform "elec- 
tric" field foi = const. According to ffTT]) . such a field corresponds to a static uniform 
supercurrent. We adopt the convention in which the supercurrent density is measured in 
units of — 2|e| (e is the electron charge), i.e., is given by K^VO. Then, for example, a current 
in the positive x direction corresponds to an "electric" field in the negative y direction, and 
if a vortex moves that way, far from any boundaries, the work done on it by the current will 
be positive. In the presence of boundaries, a nontrivial change of ^ at a boundary may give 
an additional contribution to the total work. 



III. RADIATION THEORY 

To describe tunneling, we switch to the Euclidean time, via 

T = it , 

CI4 = IQjQ , 

J' = ZJ° . 

The relations flTTl) take the form 

-id,e = f,y=B, (15) 

zdy9 = f,^ = E, (16) 
-id J = fy, = F, (17) 

where we have introduced shorthands B, E, and F, which will be much used in what follows. 
The Euclidean counterpart of the action f|T2l) is 

Se = ~KsJ dxdydr (b^ + + F^) . (18) 

Note that, as a result of the transition to the Euclidean time, the relations of all three 
components of f^^ to the derivatives of have acquired factors of i. As a consequence, on 
the instanton solution, B, E, and F will all be purely imaginary, and the action (llSp will 
be positive. This is in agreement with the a priori expectation that coupling to plasmons 
should suppress vortex tunneling (the suppression factor is e""^^) since, as it tunnels, the 
vortex has to drag the plasmon subsystem with it. 
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Maxwell equations (|T^ in these notations have the form 

- d^E - dyB = 2-Kr , (19) 
-OaF + d^B = 2njy , (20) 
d^E + dyF = , (21) 

while the Bianchi identity reads 

d^F = -d^B + dyE . 

Using these together, one can obtain independent wave equations for i?, E, and F. All that 
is left to choose, then, is a suitable form of the vortex current and the boundary conditions 
for the fields. 

We consider the theory on a rectangular strip of length L (0 < x < L) and width w 



(0 < y < w) and assume that the vortex motion is purely transverse: = 0. Then, the 
wave equations are 

dlB + V'^B = 27rd^jy , (22) 

dlE + V^E = 2nd^J^ . (23) 



Once solutions to these are obtained, the solution for F can be found from eq. (!20|) or (pTI) . 
except for the static uniform component. The latter is the static uniform "electric" field 
mentioned at the end of Sect. [Tll It corresponds to a steady supercurrent in the x direction, 
and that supercurrent can in principle have any value. Eventually, this component of F will 
be determined by the properties of the metastable state from which the system tunnels. 
The remaining (nonzero) components of the vortex current are 

jy = id,Y6{x - X)6{y - Y{t)) , (24) 
J' = t5ix-X)5iy-Yir)), (25) 

where Y{t) is the transverse position of the vortex. At a finite temperature T, Y{t) must 
be periodic in r with period (3 = 1/T. 

If the vortex could nucleate inside the strip, Y{t) would start at the upper edge, Y = w, 
move down to y = Fnuci, the nucleation point, and then back to Y = w. This would form a 
"bounce" .~ As we already noted, though, here we consider only supercurrents that are small 
enough for the vortex to have to tunnel the entire width w. Then, the relevant configuration 
is an instanton-antiinstanton (lA) pair: a vortex tunneling across the strip around time 
T = tq, plus an antivortex (or a vortex moving in the opposite direction) tunneling around 
T = Tq. a representative history of Y{t) is shown in fig. [H 

Turning to the boundary conditions (b.c), we note that, since there is no current through 
the edges of the strip, B and E satisfy, respectively, the Neumann and Dirichlet boundary 
conditions at y = and w. We can then define Fourier transforms with respect to y and r: 

B{x,y-T) = TY.e-'''"^My)BUx) , 

In 

E{x,y;r) = Tj2e-'''"^xi{y)EUx) , 

In 
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F(X) 




FIG. 1: The vortex's transverse position as a function of the Euchdean time. Ar denotes the 
instanton-antiinstanton separation, and tc the duration of an instanton. In the text, we also use 
the rescaled variable Y = Y/w, where w is the wire's width. 



where 'ipiiy) = cos{nly/w) and Xiiv) = sin^nly / w) , and the sum over / starts from Z = in 
the first case and from Z = 1 in the second; fin = '^imT are the Matsubara frequencies. In 
either case, the wave operator takes the form 

with 



The action (fT8l) becomes 

Se = -\KsWT r dx {BinBl^-n + ^^n^i.-n + F^Fi^.n) , (27) 

^ In ''^ 

where Q = 1 for / = and Q = | otherwise. 

To obtain the b.c. at x = 0, L, we need to specify how the sample connects to the outside 
world. Here, we consider the case when the leads are bulk superconductors. As a model 
of those, we use strips of the same width w and with same parameter g as the wire but of 
much larger stiffness, K'^ ^ Kg, and length L' ^ L. L' will eventually be taken to infinity. 

Consider the interface at a; = L. Denote the Fourier components of 9 there as 

Oln{.L) = 6in . 

Then, for any I and n that are not both zero, throughout the lead [x > L) 

Oinix) = ^^ne-'=-(^-^) , (28) 

where k'l^ > and is given by ([26]) with ^1 replaced by fll/gK^. Substituting (EH]) into the 
action of the lead, we obtain a contribution to the effective action of din- 

s'^ = \k',wT Qk'Jinh^n . 



2;2 



(26) 
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When we extremize the total action with respect to this term gives 9i^^n with a coefficient 
that grows at least as at large K'^. As a result, at large K'^, 6in are close to zero. The 
same applies at the other interface, at x = 0. 

We conclude that, in the case of bulk superconducting leads, both B and E satisfy the 
Dirichlet b.c. at either end: 

BiniO) = Bin{L) = , (29) 
E,„(0) = Eir,{L) = (30) 

[kin > 0). The absence of kin = from these conditions is inconsequential, since neither B 
nor E has an / = = mode (only F does). 

Note that the condition (129)) is not satisfied by the trial instanton configurations consid 



ered (for the Id case) in Ref. |28|. This explains the difference in the final results: Ref. |28 
finds that the phase slip rate remains finite at T ^ 0, while we find that it is exponential 
in 1/T. 

We now proceed to solving the wave equations for B and E, and determining F and the 
Euclidean action. 



IV. SOLUTION FOR THE PLASMON 

Since the current components flM|) and fl25l) have the same x-dependence, both B and E 
can be expressed through solutions to the equation 

-dlfi^ + klfin = dJ{x-X) (31) 

{kin > 0) with the Dirichlet b.c. fin{0) = fin{L) = 0. These solutions are readily found: 

. , 1 J cosh[kiniL - X)] sinh(fc/„a;) , x < X , , . 

•^'"^^ sinh(A;;„L) [ - cosh(fc;„X) sinh[A;,„(L - x)] , x > X , ^""^^ 

although, as we will see, many useful conclusions can be drawn even without using this 
explicit form. 

From now on, we assume that the width w is much smaller than the length L (i.e., 
w <^ L) and restrict attention to the / = (^/-independent) modes, which are the only ones 
that are potentially infrared sensitive. For tunneling paths such as the one shown in fig. [T], 
this makes our problem similar to the problem of a quantum phase slip in a genuinely Id 
wire. 

The field E has no / = modes, so it drops out of the subsequent discussion. The other 
fields are 

Bonix) = -2m{drY)nfon{x) , (33) 
2tt 

Fonix) = ^idrY)n[dJonix) + (5(x - X)] , U , (34) 

i In 

where Y is the rescaled transverse coordinate, Y = Y/w, and {drY)n is the Fourier transform 
of drY: 

{drY)n = f drdrYe'"'"^ = -l^nYn . (35) 
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The n = component of B is zero, while that of F is the static uniform component that 
should be determined from the properties (the winding number) of the metastable state. 

Let us take up the latter task first. In the wire, set Fqq = —il/KgiuT, where / needs to 
be determined. Then, in the leads, Fqq = —il/K'^wT. The total winding number is 

r^+L' I / L 2L'\ 

N{T)=tj_^^ F(x,r)rfx = -f — + — j+2vr[r(r)-Tro]. (36) 

The last term here is due to the n 7^ modes flM|) : 

iT e-*^"" Fo„(x)da; = 27rT ^ e-*^"^F„ = 27r[f (r) - TYq] • 

Note that the sum evaluates not to Y{t) but to Y{tJ without the zero mode. As a result, 
at a given time depends on the entire history of Y{t), in particular, on the value of the 
lA separation At. 

The initial and final states of tunneling correspond to points midway between the in- 
stanton and antiinstanton. These are the points at which the system enters and leaves the 
classically forbidden region. If the instanton and antiinstanton positions are sharply defined, 
i.e., Tc ^ Ar (cf. fig. [1]), the corresponding times are 

Ti = Tf- 1/2T . 

The winding numbers at these times can be computed from fl36|) and compared to those of 
ground states with uniform supercurrents. 

In particular, the winding number at t = Ti is the same as in the uniform ground state 
with supercurrent Igs given by 



Thus, the initial state belongs to the thermal ensemble built near that ground state. A 
similar relation (with a different ground-state current) applies in the final state, and we 
find, as expected, that the instanton describes tunneling between two thermal ensembles 
that differ by a 27r of the winding number. 

The duration of an individual instanton, tc in fig. [H is determined by the parameters of 
the sample. Meanwhile, as we will see, the lA separation Ar is controlled by the values of 
the temperature and current and becomes large when those are small. So, for calculating 
( 136|) . we can approximate drY in fl35l) as 

drY ^ -6{t - To) + 6{t - r^) . (38) 

Then, the last term in fl36|) becomes 

Y{n)-TYo = TAT , 

where Ar = rg — ro. In the limit L' — 00, eq. fl37|) gives 

tiK'wTAt 
I = I A 
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This relation between / and Igs can be used to learn how much action is contained in the 
I = n = mode. The main contribution comes from the leads, where Fqq = —il/K'^wT 
and Fgs^QQ = —ilgg/K'^wT. So, the I = n = term in the action (1271) . relative to the 
corresponding term in the ground state, equals 

So = \KwT J dx{\Foo\^ - |F,,,oor) = -27rJAr + 0{1/L') . (39) 

Note that this is the only term in the action that distinguishes between direct and reverse 
processes: it would change sign if we considered an antivortex tunneling at r = Tq (or a 
vortex tunneling in the opposite, positive y, direction). 

Turning to the n ^ modes, given by fl33l) and flM|) . we find that their action is 



S, = 2n^K,wT J2 [ dx {nlf^ + [djon + 5{x - X)]^} 



2 



r 

Integrating by parts in the derivative term and using eq. (!3T]) . we bring this to the form 

Using the explicit form of /on, we find for the term in the braces 

dxfon{0) + fon{x) = . , " ^, cosh[fc„(L - X)] cosh(A;„X) , 

Jo smh(/e„Lj 

where fc„ = |r2„|. This has a minimum at X = L/2 (albeit a shallow one at small /c„), 
meaning that the vortex prefers to tunnel in the middle of the wire. Setting X = L/2, we 
finally obtain 

Si = 7i^K,wTY,Kcoth{knL/2)\Y^\^ . (40) 

The approximation (1551) correctly reproduces Fourier components of drY with <^ 
I/tc- For our purposes, it will be sufficient to consider only these /c„, so we can make the 
replacement 

\Yn\' = -^[l - COs{QnM] (41) 

in fHOj) and cut off the sum in the ultraviolet at /c„ = 1/tc- 

The resulting expression for the Euclidean action, Se = Sq + Si, is applicable both to 
short (LT <^ 1) and to long {LT ^ 1) wires. In the latter case, it reproduces the result of 
Ref. Il8|, obtained by considering phase slips directly in the one-dimensional theory. In what 
follows, we restrict attention to the former case. 

We now need to determine At. We begin by considering the action PUI) for different 
values of it. For 

|Ar| < 1/T , (42) 
the sum in (HO!) can be approximated by an integral: 

/•l/rc (Iff; 

Si ^ 27rKsW / -r[^- cos(A;Ar)] coth(A;L/2) . (43) 

Jo K 
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We see that substantial contributions can only come from k ^ l/|Ar|. It is therefore 
convenient to separate the range (142|) into two. For | Ar| -C L, the cotangent can be replaced 
with unity, and we are back to the expression^*^ for the long-wire case. The underlying physics 
is that the typical wavenumber of plasmons in the final state is of order l/|Ar|, and when 
this is much larger than 1/L the effect of the boundaries is insignificant. In the present case, 
however, this is possible only for relatively large currents. Indeed, extremizing the total 
action 5*0 + 5*1 with respect to At, we obtain the saddle-point value At = Kgw/I. For this 
to be much smaller that L, we need / S> Kgw/L. Although this condition does not look 
altogether prohibitive, from now on we concentrate on the opposite, small-current, regime: 

I < . 

Then, there is no saddle point either for |Ar| ^ L or for |Ar| ~ L. 

As |Ar| increases past L — while still obeying (H2i) — the dependence of (H3l) on |Ar| be- 
comes linear, and one can verify that, for small currents, this again precludes a saddle point. 
So, we turn to At ~ 1/T and the general expression fHOl) for 5*1. 

Eq. fHOl) has an extremum at At = 1/2T (derivative of each term in the sum vanishes 
individually). This extremum is in fact a maximum, as it should be: the original integration 
was over real At = — iAr, and a saddle point that is a minimum in the real At direction is 
a maximum in real At. The contribution from 5*0 will displace the maximum from exactly 
1/2T, but for small currents this displacement is small, and for our purposes At = 1/2T 
is a good approximation. Substituting it into eq. (HTj) . we find that only odd n contribute, 
and 5*1 becomes, to logarithmic accuracy. 

Si = 87r^ KswT ■^coth{KL/2)^ —-^ + 27rKsw\n— . (44) 

n=l,3,... '^n 2L1 Tc 

where nc ~ 1/tcT. In what follows, we retain only the first, leading, term on the right-hand 
side. 

The power dissipated by vortex tunneling is given by the energy 27r/ that a vortex releases 
from the supercurrent, times the difference between the rates of the direct and reverse 
processes: 

P(J)= 27^/[7^+(/)-7^_(/)]. 

To the exponential accuracy, 7^+(/) ~ exp[— 50(1) — Si\, and = 7?.+ (— /). Expanding 

in small /, we obtain the resistance: 

RiT) ~ e-^^ ~ exp (-^) • (45) 

The activated behavior of the resistance can be interpreted by looking at the gradient 
of the phase 6 at the entry point of tunneling, t = Ti. Using the approximation fl4ip and 
setting At = 1/2T, we obtain 

dJ{x,T,)=tF{x,T,) = -^ + 4 -sm{7rn/2)[dJon{x)+5{x-X)]. (46) 

For low-frequency modes, those with < /c„ <C 1/L, 

dJon{x) + 6{x-X)^l/L, 
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so the total contribution of these modes to (H6|) is {dxO)iovi ~ vr/L. Then, the gradient energy 
contained in these modes is 



in precise accord with fH5l) . This means that the main effect suppressing resistance at low 
temperature is the population of the initial tunnehng state, due to the large gradient of 
6 already required by that time. There is, of course, an additional action associated with 
the tunneling itself, but the precise agreement between fH5|) and (jUj) implies that it is only 
subleading. 

We should note that the entirely classical appearance of the exponent in fllS]) (it does 
not require any powers of h) does not contradict it being a consequence of tunneling, rather 
than a classical, over-barrier process. Indeed, restoring h and Cq in the short-wire condition, 
under which (H5ll applies, we obtain T <^ hco/L, and this cannot be realized outside of 
quantum mechanics. 

For comparison, let us list nucleation energies for two purely classical processes. The nu- 
cleation energy of a vortex is of order ttKs (times a logarithm), and that of the independent 
saddle point, analogous to the LAMH saddle point in Id wires,— is of order Kgw/C,, where 
^ is the Ginzburg-Landau coherence length. As long as w and ^ are both much smaller than 
L, either of these energies is larger than the energy fHTl) . 

Due to a relatively large numerical factor (vr^/2) in (H7|) . the condition ^ -C L, under 
which the energy (jUj) is smaller than the LAMH activation energy, may in fact mean that 
L must be several times larger than ^. Moreover, even if this condition is satisfied at low 
temperatures, it breaks down in a region close to the critical temperature Tc. The short-wire 
condition also breaks down near in a region near Tc, since Kg and hence cq are small there; 
in that region, subleading terms in the instanton action become non- negligible. Outside of 
these regions, however, vortex tunneling is the dominant resistive process. 

For a superconductor in the dirty limit, the exponent in fHSl) can be expressed entirely in 
terms of the superconducting gap A = A(T) and the normal-state resistance Rn = pL/w, 
where p is the sheet resistivity. Indeed, in this case Kg = {nA/Ae^p) tanh(A/2T),^^ so 



where Rq = 7r/2e^ = 6.5 kQ. This suggests that the value Rn = Rq may have a special 
significance in short wires. Experimentally, it does: this is the value near which one observes 
a superconducting-insulating transition.-ii^ 



We have described a general method for calculating the effect of plasmons on vortex 
tunneling in superconducting wires and applied it to the limit of small temperatures and 
currents, when the vortex has to tunnel the entire width of the wire. The method is based 
on a duality map, through which vortices become charges and plasmons become "photons" . 
We have found that, if plasmons cannot easily leave the tunneling region, as is the case when 
the wire is short and the leads are bulk superconductors, the suppression of the resistance 
in the superconducting state is exponential at low temperatures and expressed by eq. (H5l) . 




(47) 




(48) 



V. CONCLUSION 
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As the width of the sample is made smaller and approaches the coherence length ^, vortex 
tunneling crosses over to quantum phase slips in a genuinely Id geometry. The restriction 
to modes independent of y (the transverse coordinate) that we made in sect. HVl effectively 
brings us one dimension down, so we expect the Id case to be similar to ours. 

Despite the exponential suppression, the resistance due to (thermally-assisted) vortex 
tunneling is larger, over a broad range of parameters, than that due to classical, over-barrier 
processes, such as motion of a thermally nucleated vortex-antivortex pair or a thermally- 
activated phase slip. In addition, it has a characteristic dependence on the length L of the 
sample or, equivalently, on the total normal-state resistance Rjy, cf. eq. PHl) . We hope that 
these features will allow one to distinguish between the quantum and classical processes in 
the experiment. 

The author thanks A. Bezryadin for useful comments. 
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